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Resumen 


Blanco, M. E., Vaquera, H., Villaseñor, J. A., Valdez-Lazalde, 
J. R, € Rosengaus, M. (noviembre-diciembre, 2014). 
Metodología para investigar tendencias espacio-temporales 
en eventos meteorológicos extremos: caso Durango, México. 
Tecnología y Ciencias del Agua, 5(6), 25-39. 


Se propone una metodología basada en la teoría de valores 
extremos para investigar cambios espacio temporales en 
la tendencia a largo plazo de las variables meteorológicas 
precipitación y temperatura. Se presenta un ejemplo de 
aplicación para el estado de Durango, México, analizando 
los máximos anuales de las variables temperatura máxima, 
temperatura mínima y precipitación ocurridos de 1961 a 
2000. Los resultados se presentan en mapas de tendencias 
de eventos extremos. Para el caso del estado de Durango 
se observan cambios en el comportamiento de eventos 
extremos, manifestándose una clara tendencia ala ocurrencia 
de años más secos y más cálidos en algunas zonas del estado. 


Palabras clave: series de tiempo, teoría de valores extremos, 
modelo lineal vectorial generalizado, interpolación espacial. 


Introducción 


La variabilidad temporal del clima, en particular 
de la intensidad y frecuencia de fenómenos 
climáticos extremos, puede afectar a la sociedad 
con más fuerza que los cambios en los promedios 
de dichos fenómenos (Katz € Brown, 1992). En 
el mundo, los eventos climáticos extremos son 
causa de desastres naturales (inundaciones, 
sequías, etc.), 
severamente el ambiente y la vida humana (Min, 


incendios, lo cual afecta 
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A methodology based on extreme value theory is proposed to 
investigate spatio-temporal changes in the long-term trend of the 
meteorological variables precipitation and temperature. We present 
an application for the state of Durango Mexico, analyzing maximum 
annual values of maximum temperature, minimum temperature, 
and precipitation registered during the time period 1961-2000.The 
presented results are trend maps of extreme events. For the case of 
Durango State, time changes in the behavior of extreme events are 
present, showing a clear trend towards the occurrence of dryer years 
and warmer temperatures in some areas of the state. 
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Zhang, Zwiers, £ Hegerl, 2011). En México, los 
estados de Coahuila, Chihuahua, Durango, 
Zacatecas, San Luis Potosí Tamaulipas, 
Aguascalientes, Sonora y norte de Guanajuato 
fueron declarados zonas de desastre por una 
inusual sequía durante los años 2011 y 2012. 
Según Allan y Soden (2008), y Min et al. (2011), 
es probable que la inusual frecuencia de eventos 
meteorológicos extremos esté relacionada con el 
incremento de los gases de efecto invernadero 


emitidos a la atmósfera por la actividad humana 
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y que causan el calentamiento global. La 
predicción de este tipo de eventos en el ámbito 
regional presenta un reto considerable; sin 
embargo, el esfuerzo bien vale la pena, dado 
lo esencial de la predicción en la formulación 
de estrategias para la adaptación y mitigación 
de los efectos negativos provocados por dichos 
cambios (Allan, 2011). Contar con métodos 
estadísticos adecuados para analizar el 
comportamiento espacio-temporal de variables 
meteorológicas extremas (temperatura máxima 
y mínima, precipitación máxima y mínima, 
velocidad máxima de vientos, etc.) es de 
suma importancia para la implementación de 
acciones de mitigación de dichos eventos. 
Recientemente se ha realizado un esfuerzo 
científico importante (Rahmstorf £ Coumou, 
2011) para 


temperatura inusualmente calientes (récord) 


analizar las mediciones de 
y su relación con las tasas de calentamiento 
global y de eventos de precipitación extremos 
(She et al., 2013). De acuerdo con varios autores 
(Katz, Parlange, €: Naveau, 2002; Cooley, 2009; 
Min et al., 2011; Rahmstorf € Coumou, 2011), 
el modelo probabilístico de la distribución 
de Valores Extremos Generalizada (DGEV) 
ofrece una base firme para modelar extremos 
climáticos. Una amplia base teórica sobre el 
modelo de DGEV se puede encontrar en Rolf- 
Dieter-Reiss (1997) y Coles (2001). Por otro 
lado, Yee y Stephenson (2007) proponen un 
caso particular del VGAM (Vector Generalized 
Additive Models), el VGLM (Vector Generalized 
Linear Models), para modelar los parámetros 
de la distribución de valores extremos como 
funciones lineales de covariables. 


Objetivo general 


Proponer una metodología estadística para 
identificar las tendencias espacio-temporales 
extremas delas variables climáticastemperatura 
máxima (Tmáx), precipitación máxima (P) 
y temperatura mínima (Tmín), utilizando 
conjuntamente el modelo de DVEG, el modelo 
VGAM y herramientas geoestadísticas. 


Objetivos particulares 


e Investigar si las series climáticas del estado 
de Durango presentan comportamientos 
inusuales. 

e Elaborar y mostrar mapas de tendencias 
espacio-temporales en variables climáticas 
del estado de Durango, como posible 
evidencia de áreas en donde pueden 
presentarse desastres naturales. 


Metodología 
Base de datos 


La fuente de información utilizada en este 
trabajo es la malla uniformemente distribuida 
Maya v1.0, obtenida de las interpolaciones 
de CLICOM (México Climatological Station 
Network Data) del Servicio Meteorológico 
Nacional. 

Se ubican, según coordenadas (latitud, 
longitud), los 270 nodos del estado de Durango 
de la Maya v1.0, que registran la información 
diaria de las observaciones de las variables 
meteorológicas (temperatura máxima (Imáx), 
temperatura mínima (Tmín) y precipitación 
máxima (P), del estado de Durango, en el 
periodo 1961-2000. 

Se ordena la información de cada nodo 
por fecha y se obtienen los valores extremos 
(máximos, mínimos, percentiles), por bloques 
anuales, según el teorema de Fischer-Tippet 
en la teoría de valores extremos; los valores 
máximos (mínimos, percentiles) se modelan 
con la Distribución de Valores Extremos 
Generalizada (DGEV) como sigue: 


Modelo estadístico propuesto 
Distribución de valores extremos 


Sean X,,..,X, variables aleatorias indepen- 
dientes e idénticamente distribuidas, con 
función de distribución F. Sea el máximo 
de las observaciones M, = máxiX,,...,X ). El 


teorema de Fischer Tippet (Rolf-Dieter Reiss, 
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lim 


n—0w 


1997) indica que 


M4, 
pita <=)-< (>) 


n 


con (a, |ER y [b, >0); Gíz) es llamada la 
Distribución de Valor Extremo Generalizada 
(DVEG), definida por: 


cap (eo 


G(z)=exp exp al ,Ex0,1 


al 


N 
I 


+E2)>0 


10] 


(1) 


Donde u es el parámetro de localidad, o es 
el parámetro de escala y E es el parámetro de 
forma. 

Según el límite de DGEV (1), cuando E = 0, 
tenemos la llamada distribución Gumbel, la 
distribución Fréchet con E > 0 y la distribución 
Weilbull con 8 < 0. 


Modelo de regresión usando la DGEV 


Un modelo que puede usarse para realizar 
regresión en los parámetros de DGEV (1) 
es el Modelo Lineal Generalizado Vectorial 
(VGLM), un caso particular del Modelo Lineal 
Generalizado de Vectores Aditivos (VGAM) 
propuesto por Yee y Stephenson (2007). El 
modelo VGLM permite que los parámetros de 
la distribución de valores extremos se puedan 
modelar como funciones lineales o suavizadas 
de covariables. Esto es: 
y¡(x, —DGEV ((x,),0(x;),8(x,)) (2) 
La regresión en (2) de los parámetros del 
modelo DGEV (1) se realiza como sigue: 


0 +Iog[E+3 pr +Bos 6) 


n= log (0) =P, +B2X; (4) 


n¿=u4=By (5) 


con coeficientes de regresión desconocidos y 


donde  y= (mast5,1,) = (tos(0).los(E+ ) 


de funciones 
monótonas, llamadas “función liga”. Es posible 


representa la coordenada 
estimar los parámetros de la distribución 
DVEG y del modelo VGLM usando un paquete 
computacional del programa R, 
VGAM (Yee, 2008). 

Se utiliza una prueba de razón de vero- 
similitudes 


llamado 


generalizada para probar la 
estacionariedad, usando la verosimilitud de los 
modelos completo y reducido. 


Prueba de razón de verosimilitudes 
generalizada para tendencia en la serie de 
extremos 


Se investiga la existencia de tendencia en series 
de valores extremos, se incluye el factor tiempo 
(£) como covariable x, en la regresión (3), (4), 
(5) de los parámetros del modelo DGEV (1). 
Así, para Y, - DGEV(wolt), E(t)), los valores 
máximos (mínimos) anuales tomados en el 
tiempo (t). Se plantea la hipótesis de interés, 
como sigue: 

H,; no existe evidencia estadística de 
tendencia en los parámetros de la distribución 
de valores extremos generalizada: 


Modelo reducido (My) 

= Y, —GEV(n3"(B)m2*(B2),0m5*(8)) 0 

H,: existe evidencia estadística de tendencia 
con respecto a los parámetros de escala y de 
forma: 


(7) 
= Y, —GEV(n,* (Bs1).17* (Ba, + Bast),n;* (Br + Brat) 


Modelo completo (M, ) 


Se utiliza la estadística de prueba: 


Xo =2|1(M,)-e(Mp)]-x%-», (8) 
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Donde p,es el número de parámetros de M; 
1(M,) es la log-verosimilitud del modelo M,con 
¡=1,0. 

La regla de decisión es rechazar la hipótesis 
nula si xj es mayor que un valor de tablas de la 
distribución x;,,yy Se utiliza el paquete VGLM 
bajo el programa R (Team, 2008) para ajustar 
el modelo completo (7) (cuando el tiempo 
es la variable predictora), donde se asume 
que la serie es no estacionaria y suponiendo 
estacionariedad en el modelo reducido (sin el 
tiempo como covariable). Los coeficientes f,, 
y f,, del modelo (7) representan el efecto de 
tendencia en la variable de extremos. 

Se estiman los coeficientes de tendencias 
bajo el modelo VGLM, los cuales, dado que 
se obtienen usando el método de máxima 
de verosimilitud para muestras grandes, se 
distribuyen normalmente. 


Método de Interpolación 


Se aplica el método kriging ordinario para 
interpolar las tendencias. Se presentan los 
mapas de tendencia para las variables meteo- 
rológicas utilizando la aplicación ArcGIS 8.1. 

Kriging es un método de interpolación 
espacial usado en investigación climática 
(Drignei, 2009). El modelo ordinario de kriging 
es el siguiente: 


Z(s)=u+e(s) (9) 


donde s = (x,y) es un punto en el espacio. El 
predictor está formado por la suma ponderada 
de los datos: 


N 
Z (sp) = Y r,Z(s;) (10) 
¡A 
donde s, es el punto sobre el que se requiere la 
estimación. N es el número de puntos alrededor 
de s, que se usan para la predicción. Los A, son 
los pesos para cada una de las mediciones. 
Z(S,) son los valores observados en el punto S, 
Z(S,) es la estimación sobre S,, 
Los estimadores A, son los valores que 
minimizan: 


Z(sp)= Saz] (11) 


Sujetos a la restricción de que la suma de los 
pesos A, sea igual a 1. 

La solución a la minimización, sujeto a la 
restricción, está dada por la ecuación kriging 
siguiente: 


A=T*g (12) 
Donde: 
Yu Yin 1 A, Yio 
Pe E a eS 
Ya + Yu 1 Ay Ywo 
1..1.0 m 1 


Donde m es desconocido y es estimado 
de forma que se cumpla con la restricción 


de > + Y 


ds con a semivariograma ajustado, 


son las semivarianzas 


con la distancia entre los puntos s, y s; y los 
Y, Son las semivarianzas calculadas con el 
semivariograma ajustado, con la distancias 
entre el punto s, y los puntos s.. 

El semivariograma ajustado es el modelo 
que ajusta al semivariograma empírico, el 
cual es una gráfica entre las distancias y la 
semivarianza promedio observada en dicha 
distancia. La semivarianza entre dos puntos es 
igual a 0.5 multiplicado por la diferencia entre 
los valores elevada al cuadrado. 

El modelo esférico para el semivariograma 
es definido como sigue: 


AECI 
1(1,0)=4 *[2 0, (28, 


O, para 0, < 11d (13) 


para O< [la <0, 


Donde 6, > 0 es el parámetro de umbral y 
9, > 0 es el parámetro de rango. 
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Validación del modelo 


La validación del modelo es usando los valores 
de log verosimilitud reportados por VGAM. 


Descripción esquemática 
En las figuras 15 y 16 se muestra la síntesis de 


la metodología propuesta y el diagrama en 
general de la propuesta para analizar datos de 


Resultados 


Análisis descriptivo de las variables 
temperatura máxima (Imáx) temperatura 
mínima (Imín) y precipitación máxima (P) 


La variable Tmáx muestra una marcada 
variación en las pendientes de regresión 
asociadas con el parámetro de escala, lo cual 
indica que a lo largo de los cuarenta años ana- 
lizados, los valores se han alejado de la media 


eventos extremos, con la finalidad de que el en temperatura máxima, significativamente. 


La ubicación de los 270 nodos (Maya v1.0) y de 
algunas localidades (puntos grises) en el estado 


uso de esta metodología aporte información, 
con otro conjunto de datos, en el análisis de 


valores extremos. de Durango se muestran en la figura 1. 


Maya v1.0, con 270 nodos para Durango, México, que contienen 
series de tiempo de observaciones diarias, con comportamientos 
inusuales 


Base de datos Se ordenan los 270 nodos por fecha y se obtienen bloques anuales 


Se extrae de los bloques anuales los valores máximos, (Tmáx y P), 
mínimos (Tmín), percentil 10 (Tmín) y percentil 90 (Tmáx y P) 


Los valores máximos, mínimos y percentil (Tmáx, Tmín) por 
bloques anuales se modelan con la DGEV (1) 


Para identificar tendencias en los parámetros de regresión la 
DGEV (1) se propone el VGLM 


Modelo 
propuesto 
Metodología 
propuesta Con VGLM se estiman los parámetros de regresión en los 
parámetros de la DGEV (1), en función del tiempo modelo completo 
(7) y modelo reducido (6) de la prueba de razón de verosimilitudes = 
generalizada S 
3 
2 
a 
E 
3 
ñ Se obtienen mapas de tendencias con los coeficientes de regresión 07 
: Método de ma ASOCIadOs con los parámetros de escala y forma de la DGEV, con el E 
interpolación método Kriging S 
E 
a 
Materiales ma Paquete computacional VGLM del programa R y la aplicación > 
ArcGIS 8.1. de 
2 
3 
Figura 15. Descripción esquemática de la metodología que se propone. - 
3 
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* Bloques anuales de máximos, 
mínimos y percentiles (se distribuyen 


DGEV) 


suas e Ubicación de 


posibles zonas 
propensas a 
desastres naturales 


e Una alternativa 
para identificar 
tendencia en 
eventos extremos 


Regresión en los parámetros de la 


e Eventos extremos en 
comportamiento de series 
de tiempo de variables 
meteorológicas 


e ¿Tendencia de eventos 
extremos? 


DGEV 


Metodología propuesta 


DGEV con el modelo propuesto 
VGLM con Tmáx, Tmín y P como 
variables respuestas 


Mapas de tendencias 


Prueba de razón de verosimilitud para 
obtener coeficientes de regresión 
asociados con los estimadores de los 
parámetros de forma y escala de la 


Interpolación de coeficientes de 
regresión (tendencias) 


Figura 16. Diagrama en general de la propuesta para analizar datos de eventos extremos. 


En las figuras 2, 3 y 4 se plasma el comporta- 
miento de las series de tiempo de las variables 
climáticas temperatura máxima anual, tempe- 
ratura mínima anual y precipitación, en algu- 
nas localidades señaladas. 

Algunas localidades al norte del estado 
de Durango presentan temperatura máxima 
mayor que algunas localidades al sur; El Salto 
registra la temperatura máxima menor. 

En la figura 3, algunas localidades al sur 
del estado de Durango presentan temperaturas 
mínimas menores que algunas localidades al 
norte a lo largo de cuarenta años; la localidad 
El Salto registra las temperaturas mínimas 
menores. 

En la figura 4 se indica la variación anual 
en la precipitación máxima (milímetros) de 
las localidades. Las localidades de El Salto y 
Durango presentan en el incio de la tercera 


década una tendencia positiva; sin embargo, es 
complicado observar alguna tendencia. 

En la figura 5 se presenta el comporta- 
miento de los valores máximos y mínimos de 
temperatura (promedio diario de 14 600 días), 
del percentil 90 (Tmáx, P) y del percentil 10 
(Tmín). Los valores representan el promedio 
diario de los 270 nodos. 

Los valores de percentil y los valores de 
máximos de los promedios diarios de las 
variables (Imáx, Tmín, P) indican que no 
hay tendencia, aunque el comportamiento 
de la serie de tiempo en las observaciones 
de percentil muestra que es posible cometer 
errores; a pesar de los cambios bruscos que se 
observan, ello no es evidencia de tendencia 
con los valores máximos. Luego es preciso 
realizar un estudio por nodo, en vez de los 
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Figura 1. Nodos de Maya v1.0 ubicados en el área de estudio. 
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Figura 2. Series de tiempo anuales de Tmáx (grados) de algunas localidades del estado de Durango (1961-2000). 
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Figura 3. Series de tiempo anuales de Tmín (grados) de algunas localidades del estado de Durango (1961-2000). 
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Figura 4. Series de tiempo de precipitaciones máximas (milímetros) anuales de algunas localidades del estado 
de Durango (1961-2000). 
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Figura 5. Series de tiempo de valores máximos y percentiles de promedios diarios, anuales de 270 nodos del estado 
de Durango (1961-2000). 
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Figura 6. Series de tiempo de valores máximos de Tmín registrados en las bases de datos CLICOM y MAYA v1.0 
para Cuecamé y El Salto, Durango, periodo (1961-2000). 
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Figura 7. Series de tiempo de valores máximos de precipitación máxima registrados en las bases de datos CLICOM 


y Maya v1.0 para Cuencamé y El Salto, Durango, periodo 1961-2000. 
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valores promedios diarios como gráfica repre- 
sentativa del estado de Durango. 

Las figuras 6, 7 y 8 presentan un comparativo 
de las series de tiempo registradas en las bases 
de datos CLICOM y Maya v1.0, para algunas 
localidades. Se usaron valores máximos. 

En la figura 6, la serie de tiempo registrada 
en Maya v1.0 para la localidad de Cuencamé 
presenta similitud en su comportamiento, 
comparada con la serie de tiempo de CLICOM. 
Las series registradas en CLICOM para la 
localidad El Salto no está completa. 

En la figura 7 se observa la serie de tiempo 
de la localidad Cuencamé; presenta valores 
máximos de Maya v1.0 por debajo de las 
observaciones de CLICOM, con diferencias 
significativas. El Salto carece de información al 
término de la tercera década. 

En la figura 8 se indica la serie de tiempo de 
la localidad de Cuencamé, que registra valores 
mayores de Maya v1.0 en relación con los de 
CLICOM; los valores máximos de El Salto de 
CLICOM por encima de los valores máximos 
de Maya v1.0. 

De las figuras 6, 7 y 8 se nota que es posible 
tomar un valor extremo en Maya v1.0, que esté 
muy alejado del valor extremo de CLICOM, para 
ubicar si los posibles cambios en las tendencias 
se preservan; se comparan los valores máximos 
anuales de las variables meteorológicas (Tmáx, 
Tmín, P) con los valores percentiles 90 anuales 
(en Tmáx, P) y percentil 10 anual (en Tmín), y se 
implementa la metodología propuesta tanto en 
valores máximos anuales como en percentiles 
anuales de Maya v1.0. 


Estimación de tendencias en el 
comportamiento de las variables climáticas 


Los cuadros 1, 2 y 3 contienen el modelo 
completo (7), el modelo reducido (8) y las 
pendientes de regresión (3) y (4) asociadas con 
los parámetros de escala y forma del modelo 
(1), además de los resultados de la estadística 
de prueba (8). 

En el cuadro 4 se presentan las tendencias en 
el ajuste de VGLM con respecto a los parámetros 


de escala y forma en los 270 nodos en el estado 
de Durango, bajo la hipótesis nula (tendencia) 
y la hipótesis alternativa (no tendencia) de la 
prueba de razón de verosimilitudes. 

En el cuadro 5 se muestran las tendencias 
en el ajuste de VGLM, en cuanto a los valores 
percentil 90 y 10 de los coeficientes de regresión 
asociados con los parámetros de escala y forma 
en los 270 nodos en el estado de Durango, bajo 
la hipótesis nula (tendencia) y la hipótesis 
alternativa (no tendencia) de la prueba de 
razón de verosimilitudes generalizada. 


Mapas de tendencias en las variables 
climáticas. 


Se usan los coeficientes de regresión P,, y 
B,, (las pendientes de regresión asociadas 
con los parámetros de escala y forma, 
respectivamente), para obtener los mapas de 
tendencia para las variables meteorológicas 
mediante interpolación en la aplicación ArcGIS 
8.1. Los resultados se presentan a continuación 
(figuras 9 a 14). 

Se pudo notar que hay similitudes en ciertas 
regiones al comparar los mapas de tendencias, 
obtenidas con valores máximos anuales y los 
valores percentil 90 anuales. 

En las figuras 9 a 14, se observa que hay 
un aumento considerable en la tendencia 
(positiva y negativa) en el parámetro de escala, 
en temperaturas máximas anuales en algunas 
áreas del estado, en temperaturas mínimas y 
precipitación máxima anual, lo cual significa 
que estas regiones son más propensas a 
manifestar algún fenómeno climático extremo 
(por la tendencia de los cuarenta años de 
estudio). Los mapas generados se pueden 
emplear para localizar áreas donde eventos 
extremos pueden ser aún más severos en 
el futuro. Al norte y al oeste del estado se 
presentan las temperaturas máximas mayores; 
al sureste del estado, temperaturas mínimas 
menores; y al sur y suroeste, precipitación 
máxima menor. 
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Cuadro 1. Modelo completo (M,), donde f,, y P,, corresponden a las pendientes de regresión asociadas con los parámetros 
de la DGEV (1), y “(M.) el valor de log-verosimililutd. 


Modelo completo (M,) 


Localidad Br Bas UM) 
El Salto 0.10045735 0.002131688 -156.5825 
Cuencamé -0.20115133 -0.006367134 -138.919 


0.30995383 


0.42849085 -149.7359 


Durango 


Cuadro 2. Modelo reducido (6) de la prueba de razón de verosimilitud (8) para probar tendencias en la precipitación 
máxima (1961-2000) en algunas localidades de Durango. 


Modelo reducido (M,) 


Localidad (MM) 
El Salto -156.9733 
Cuencamé -141.5256 
Durango -161.2376 


Cuadro 3. Prueba de razón de verosimilitud (8) para detectar tendencias en la precipitación máxima (1961-2000) en algunas 
localidades de Durango. El estadístico de prueba es xj y; = 3.841459. 
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C 


« vj30/0u32], 


Localidad Regla de decisión Tendencias 
El Salto 0.78152139 < Xdo951) 0 

Cuencamé AMES sn il 
Durango 23.003907 > Xio9s,1) 1 


Cuadro 4. Resultados de ajuste usando VGAM, con valores máximos (para Tmáx, P) y mínimos (Tmín). 


Variable Nodos no convergentes Nodos convergentes Tendencias significativas al 5% 
Precipitación máxima 20 250 93 
Temperatura máxima 18 202 161 
Temperatura mínima 11 259 157 


Cuadro 5. Resultados de ajuste usando VGAM, utilizando valores percentil 90 (Tmáx, P) y percentil 10 (Tmín). 


Variable Nodos no convergentes Nodos convergentes Tendencias significativas al 5% 
Precipitación 15 255 97 
Tmín 20 250 183 
Tmáx 12 258 196 
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Análisis de tendencia de la pendiente de regresión, asociada al parámetro Análisis de tendencia de la pendiente de regresión, asociada al parámetro 
de escala de la DGEV, Tmáx, utilizando valores máximos en bloques anuales de forma de la DGEV, Tmáx, utilizando valores máximos en bloques anuales 


Mapa de predicción Mapa de predicción 
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Figura 9. Mapas de tendencia de la pendiente de regresión asociada con el parámetro de escala y de forma de la DVEG 
para la variable Tmáx. Se usaron valores máximos. 


Análisis de tendencia de la pendiente de regresión, asociada al parámetro Análisis de tendencia de la pendiente de regresión, asociada al parámetro 
de escala de la DGEV, Tmáx, valores percentil 90 en bloques anuales de forma de la DGEV, Tmáx, valores percentil 90 en bloques anuales 


Mapa de predicción Mapa de predicción 
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Figura 10. Mapas de tendencia de la pendiente de regresión asociada con el parámetro de escala y de forma de la DVEG 
para la variable Tmáx. Se usaron valores del percentil 90. 


Análisis de tendencia de la pendiente de regresión, asociada al parámetro Análisis de tendencia de la pendiente de regresión, asociada al parámetro 
de escala de la DGEV, Tmín, valores mínimos en bloques anuales de forma de la DGEV, Tmín, valores mínimos en bloques anuales 


Mapa de predicción 
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Figura 11. Mapas de tendencia de la pendiente de regresión asociada con el parámetro de escala y de forma de la DVEG 


para la variable Tmín. Se usaron valores máximos de los datos de Maya v1.0. 
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Análisis de tendencia de la pendiente de regresión, asociada al parámetro 
de escala de la DGEV, Tmín, utilizando percentil 10, en bloques anuales 


Mapa de predicción 
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Análisis de tendencia de la pendiente de regresión, asociada al parámetro 
de forma de la DGEV, Tmín, utilizando percentil 10 en bloques anuales 
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Figura 12. Mapas de tendencia de la pendiente de regresión asociada con el parámetro de escala y de forma de la DVEG 


para la variable Tmín. Se usaron valores percentil 10 de los datos de Maya v1.0. 


Análisis de tendencia de la pendiente de regresión, asociada al parámetro de 
escala de la DGEV, precipitación, utilizando valores máximos en bloques anuales 
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-0.049252 - -0.037035 
22 -0.037035 - -0.027958 
HA -0.027958 - -0.021214 
MA -0.021214 - -0.016204 
HA -0.016204 - -0.012481 
-0.012481 - -0.007471 
MA -0.007471 - -0.000727 
MX -0.000727 - 0.008350 
0.008350 - 0.020567 
MA 0.020567 - 0.037011 


we 
S 


uu Kilómetros 
0 60120 240 


Análisis de tendencia de la pendiente de regresión, asociada al parámetro de 
forma de la DGEV, precipitación, utilizando valores máximos en bloques anuales 
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Figura 13. Mapas de tendencia de la pendiente de regresión asociada con el parámetro de escala y de forma de la DVEG 


para la variable precipitación máxima anual. 


Análisis de tendencia de la pendiente de regresión, asociada al parámetro 
de escala de la DGEV, precipitación, valores percentil 90, en bloques anuales 
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Análisis de tendencia de la pendiente de regresión, asociada al parámetro 
de escala de la DGEV, precipitación, valores percentil 90, en bloques anuales 
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Figura 14. Mapas de tendencia de la pendiente de regresión asociada con el parámetro de escala y de forma de la DVEG 


para la variable precipitación. Se usaron valores del percentil 90 de los datos de Maya v1.0. 
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Conclusiones 


Se muestra y prueba una metodología para 
investigar el comportamiento espacio-tiempo 
de valores extremos en series meteorológicas. 
Las tendencias obtenidas para el estado de 
Durango en las variables estudiadas sirven 
para monitorear los posibles cambios en 
variables climáticas y anticipar escenarios ante 
un eventual cambio climático. 
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